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Abstract 

We study the Cosmic Microwave Background (CMB) anisotropies produced by 
cosine-type quintessence models. In our analysis, effects of the adiabatic and isocur- 
vature fluctuations are both taken into account. For purely adiabatic fluctuations 
with scale invariant spectrum, we obtain a stringent constraint on the model param- 
eters using the CMB data from COBE, BOOMERanG and MAXIMA. Furthermore, 
it is shown that isocurvature fluctuations have significant effects on the CMB angular 
power spectrum at low multipoles in some parameter space, which may be detectable 
in future satellite experiments. Such a signal may be used to test the cosine-type 
quintessence models. 



1 Introduction 



It is widely believed that the present cosmological observations are consistent with low 
density Cold Dark Matter (CDM) models with scale invariant adiabatic density fluctu- 
ations. The required energy density of non-relativistic matter is about 30 — 40% of the 
critical density (i.e., fi m ~ 0.3 - 0.4). On the other hand, the recent BOOMERanG 
and MAXIMA [[| experiments on anisotropies of Cosmic Microwave Background (CMB) 
strongly suggest that our universe is flat, which is also the prediction of the inflationary 
universe. Thus, there exists dark energy which fills the gap between Q m and the total 
energy density O tot . Although usually the dark energy is assumed to be the cosmological 
constant, a slowly evolving scalar field with positive energy can also account for the dark 
energy || [|. Such a scalar field is called quintessence and has been studied by many 
authors JJ. 

At present, the quintessence models are classified into two types; tracker type || and 
cosine type ||. The former has an attractor-like solution which explains the present dark 
energy without fine-tuning of the initial condition, while the latter needs tuning of the 
initial value of the scalar field. This is why the tracker type is favored among cosmologists. 
However, even for the tracker-type models the model parameters should be fine-tuned to 
produce the required value of the present quintessence density. In addition, the potential 
for the tracker field is very exotic and hard to realize in particle physics models.^ 1 ] 

One of the observational effects produced by the existence of the quintessence is the 
CMB anisotropies. In many cases, the quintessence dominates the universe at late times 
(z ~ 0(1), with z being redshift) after the recombination. At that epoch the gravitational 
potential is changed because the equation of state for quintessence is different from that 
for non-relativistic matter, which leads to an enhancement of the CMB anisotropies at 
large angular scales / ^ 10 (where / is the index of spherical harmonics) due to the late- 
time integrated Sachs- Wolfe effect ]7|. The quintessence also changes the locations of 
the acoustic peaks in the CMB angular power spectrum || because the projection of the 
horizon at last scattering onto the present sky is enlarged compared with models with the 
cosmological constant. 

Furthermore, in the case of cosine-type models, during inflation the quintessence scalar 
field is effectively massless and has significant quantum fluctuations as large as H- in t/2ix 
where ifj n f is the Hubble parameter during inflation. These fluctuations behave as isocur- 
vature mode and hence the quintessence may have both adiabatic and isocurvature per- 
turbations. The isocurvature fluctuations produce adiabatic ones once their wavelengths 
enter the horizon. Therefore, we expect larger anisotropies in CMB angular spectrum for 
cosine-type quintessence mo dels, p 7 

In this paper we study the CMB anisotropies produced by quintessence models paying 
attention to the cosine-type quintessence models. It is shown that the isocurvature fluc- 

# 1 For model building of the tracker field, see Refs. ||. 

# 2 Abramo and Finelli || studied the isocurvature fluctuations for tracker-type quintessence models 
without considering the evolution of the tracker field and its fluctuations during inflation. 
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tuations have significant effects on the CMB angular power spectrum at low multipoles in 
some parameter space, which may be detectable in future satellite experiments and can be 
an interesting signal of the cosine-type quintessence models. 

This paper is organized as follows. In Section 0, we discuss dynamics of the scalar-field 
zero mode. In Section behaviors of the fluctuation in the quintessence amplitude are 
discussed. Then, in Section [|, CMB anisotropy in cosmological models with quintessence 
is numerically calculated. Section || is devoted for conclusions and discussion. 

2 Dynamics of Scalar-Field Zero Mode 

We start our discussion with the behavior of the zero mode of the quintessence field Q(t, x). 
We denote the zero mode of the quintessence field as Q(t). 

The energy-momentum tensor of the quintessence field is given by 

= d^QdvQ - ~ [d a Qd a Q + 2V(Q)] g^, (2.1) 

where V is the quintessence potential. Thus, for Q = Q, the energy density and the 
pressure of the quintessence field are 

Pq = Iq+V(Q), Pq = \q-V{Q) 1 (2.2) 

respectively, where the "dot" represents the derivative with respect to time t. If the kinetic 
energy of Q is negligible compared to the potential energy, the equation of state becomes 
ujq = pq/pq ~ —1 and the energy-momentum tensor of the quintessence field behaves like 
that of the cosmological constant. 

The zero mode Q obeys the following equation of motion 

Q + 3HQ + V\Q) = 0, (2.3) 

where the "prime" denotes the derivative with respect to Q. In the very early universe 
where H 2 ^> |V/Q|, the slow-roll condition is satisfied and the motion of Q is negligible 
even if Q is displaced from the minimum of the potential. If the slow-roll condition is 
satisfied until today, the equation of state ujq is always very close to —1 and the quintessence 
is indistinguishable from the cosmological constant. However, this is not always the case. If 
the effective mass of the quintessence may become larger than the expansion rate, Q starts 
to move and uoq varies. In addition, since the quintessence is a dynamical field, its energy 
density may fluctuate. Due to these facts, quintessence models may have an interesting 
consequence in the evolution of the cosmological perturbations, in particular, in the CMB 
anisotropy. The epoch when the quintessence starts to move is strongly model-dependent. 
If the slow-roll condition is satisfied until very recently, the quintessence field may provide 
significant amount of the energy density with the equation of state close to —1 to the 
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total energy density of the universe. This may be a solution to the "dark energy" problem 
alternative to the cosmological constant. 

In this paper, we consider the cosine-type quintessence potential: 



V(Q)=A i 




^S- (2 - 4) 



This type of potential can be generated if the quintessence field is a pseudo Nambu- 
Goldstone boson. In this class of models, effective mass of the quintessence field is always 
of 0(A 2 //q), and is insensitive to the amplitude Q. Requiring that the slow-roll condition 
for the quintessence field be satisfied until very recently to realize ujq ~ —1, the com- 
bination A 2 / fq cannot be much larger than the present expansion rate of the universe. 
Consequently, motion of the quintessence field is negligible for z ^> 1. Thus, in this class 
of models, the present energy density of the quintessence field sensitively depends on the 
initial amplitude. 

In the cosine-type models, the present density parameter for the quintessence depends 
on the following three parameters: Jq, A, and the initial amplitude of the quintessence 
field denoted as Q\. Assuming the flat universe (i.e., Q m + Qq = 1), however, one relation 
among these parameters holds once we fix the present matter density Q m = Qcdm + Qb, 
where ^cdm and Q\, are density parameters for the CDM and baryons, respectively. 

In Fig. [I], we plot the contours of constant Q\ which realizes Qq = 0.7 in the flat 
universe. Here, we take h = 0.65, where h is the present Hubble parameter in units of 100 
km/sec/Mpc. Notice that there is no possible value of Qi consistent with Qq = 0.7 for 
A < 1.9 x 10~ 12 GeV. This is because the energy density of the quintessence is at most 2A 
and hence pg cannot be large enough if 2A 4 < Qqp c with p c being the critical density of 
the present universe. When 2 A 4 is close to QqPc, Qi should be close to tt/q to realize the 
relevant value of Qq. In this case, motion of the quintessence field is almost negligible and 
energy density of the quintessence behaves like that of the cosmological constant. 

On the contrary, for larger A, Qi smaller than ttJq is possible and the motion of the 
quintessence becomes important. In particular, if A becomes large enough, the quintessence 
oscillates around the minimum of the potential. As one can see, for large enough A, the 
contours have an oscillatory behavior. This can be understood as follows. When the 
relation H ^ A 2 / Jq is satisfied, the quintessence field starts to oscillate. If the combination 
A 2 //q is large, the oscillation starts earlier epoch and the quintessence field undergoes many 
oscillations until the present time. We can read off this behavior from Fig. |I[ We call 
the parameter space where the oscillation of the quintessence is significant as "oscillatory 
region." When the quintessence field starts to oscillate earlier, it dominates the energy 
density of the universe from earlier epoch. This has significant implications to the CMB 
power spectrum. 

Evolution of ujq is also important since it can be quite different from the case of the 
cosmological constant, in particular when the motion of the quintessence is non- negligible. 
Typical behaviors of ujq as a function of the scale factor are shown in Fig. |2|. When the 
curvature of the potential is smaller than the expansion rate of the universe, evolution 
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Figure 1: Contours of constant Qi/ Jq = 3.0,2.75,2.5,2.25,2.0,1.75 and 1.5 from left to 
right. The value of Qi is determined such that Qq = 0.7 in the flat universe (i.e., fl m + 
Qq = 1) today. The Hubble parameter is taken to be h = 0.65. The parameter space 
A<1.9x 10~ 12 GeV is the region where f2q = 0.7 cannot be realized for any values of Q\. 



of the quintessence can be neglected, and ujq ~ —1. On the contrary, when H^zA 2 /fq, 
ujq oscillates between —1 and +1. If the oscillation is fast enough, ojq is effectively 
taking the average over the oscillation. In this case, the scalar field Q behaves as a non- 
relativistic matter. If the condition H ~ A 2 //q is realized in a very recent universe, 
however, the situation is quite different. In this case, the quintessence field starts to move 
very recently and the averaged equation of state can be smaller than 0. If so, the scalar field 
Q behaves differently from non-relativistic matter and cosmological constant. This may 
have interesting impacts on the CMB anisotropy as we will see in the following sections. 
We also present the evolution of the density parameter of the quintessence: 

P Q (t) 



PtotW 



(2.5) 



Like the case of the cosmological constant, fig(t) is negligibly small in the early uni- 
verse. When z ~ 0(1), however, fig(t) becomes close to 1 and the energy density of the 
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Figure 2: Evolution of the equation of state ujq for (a) /q = 5.0 x 10 17 GeV and A = 
2.3 x 10" 12 GeV (solid line), and (b) f Q = 5.0 x 10 17 GeV and A = 4.0 x lO" 12 GeV 
(dashed line). Here, we take h = 0.65, f2 m = 0.3 and Qq = 0.7. 



quintessence becomes sizable. One interesting feature is that, if the slow-roll condition does 
not hold, the behavior of ^q(£) is quite different from that of the cosmological constant. In 
particular, as A becomes larger, pq becomes the dominant component of the total energy 
density at earlier stage. 

The observations of Type la Supernovae suggest that ojq is less than about —0.7 at 
z^S 1 for fi m = 0.3 [TIJ. However, it was pointed that there exist some possible systematic 
errors such as dust absorption and/or evolution effect [11|. Therefore, in this paper, we do 



not use the supernovae data to obtain constraints on the quintessence models. 



3 Fluctuation in the Quintessence Field 

Since the quintessence is a scalar field, its amplitude may have position-dependent fluctu- 
ations. To investigate its behavior, we decompose the Q field as 

Q(t,x) =Q(t) +q(t,x), (3.1) 

where q is the perturbation of the amplitude of the quintessence field. Hereafter, we study 
the evolution of q using the linearized equations for the perturbations. 
The equation of motion for q is, in the synchronous gauge, 

q + 3Hq - (±) ~' d\q + V"(Q)q = ~hQ, (3.2) 
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Figure 3: Evolution of the density parameters. We show VLq (solid line), fl m (dashed line) 
and f2 radiation (dotted line). Cosmological and model parameters are the same as Fig. |2| 



where d% is the derivative with respect to the comoving coordinate x \ Here, the perturbed 
line element in the synchronous gauge is given by 



ds 2 



a 



-dt + — + hij)dx l dx 3 = — -dr 2 + (5^ + h ij )dx % dx 3 



«o 



(3.3) 



where r is the conformal time coordinate, a the scale factor at time t, cto the scale factor 
at the present time, and h is the trace of hg^ 3 In the momentum space, we expand as 



hij (t, x) 



d 3 k 
(2tt) 3 / 2 



ikx 



(3.4) 



with k 2 = kik 1 . Notice that, in the momentum space, the gauge-invariant variable $ is 
related to h and r] as 



*(*) = ^ 



2k 2 



d 2 h(k) +Q d 2 7](k) + Ida (dh(k) + Q dr](k)\ 



Or 2 



Or 2 



adr \ dr 



dr J 



(3.5) 



Numerically, we checked that the terms proportional to the derivatives of r] are subdominant 
relative to contributions from the derivatives of h. 

There are two origins of non-vanishing q at the present time. As indicated in Eq. 
( p.2|) , the metric perturbation h generates non- vanishing q even if q initially vanishes. We 



# 3 The metric perturbation h should not be confused with the Hubble parameter. 
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call this fluctuation as "adiabatic fluctuation" of the quintessence field. Let us briefly 
discuss how the adiabatic fluctuation of q is generated. For this purpose, it is instructive 
to consider the matter-dominated universe even though, in our case, the energy density 
of the present universe has significant contribution from the quintessence. Fluctuation for 
the scale smaller than the horizon scale oscillates and damps, so we consider q(k) which 
has larger physical momentum than the expansion rate of the universe. Then, the third 
term in the left-hand side of Eq. (|3.2| ), i.e., dfq, is negligible. In solving Eq. (|3.2|) , it is 
convenient to use the fact that \I> is (almost) constant of time in matter-dominated epoch. 
In addition, for the zero mode, we use the slow-roll condition to derive Q ~ — |(V/if). 
With these informations, q(k) at a given time t is, in the synchronous gauge, 

Q(k) = ^(a /a) 2 k 2 V'*(k)t 4 , (3.6) 

which is proportional to t 8 ' 3 . In the realistic situation, however, some of the approximations 
used to derive the above relation may fail. Therefore, in our analysis, we numerically solve 
the equation of motion for the quintessence field to obtain q(k). 

The second origin is the primordial perturbation in the quintessence amplitude gener- 
ated in the very early universe, probably during the inflation. We call this fluctuation as 
"isocurvature fluctuation," since the total density fluctuation and the potential ^ vanish 
as a — > if this is the only source of the fluctuation in the early universe. 

The isocurvature mode may arise due to the quantum fluctuation during the inflation. 
In order to calculate the expected fluctuation generated during the inflation, we quantize 
the scalar field in the de Sitter background. Identifying q as a field operator, it can be 
expanded as 

q(t, x) = f \a k Mr)e rk * + aM^e'^} , (3.7) 



(2tt) 3 / 2 

where, assuming minimally coupled Lagrangian for the quintessence field, the mode func- 
tion is given by 

(pk (r) = ^H ml r^Hl 1 \kr), (3.8) 

with being the Hankel function of the first kind and 

9 9 ml 

Here, m q is an effective mass for the quintessence field during the inflation, and H in { is the 
expansion rate of the de Sitter background. 

Adopting the Bunch-Davis vacuum jl2]], the operators and a k are identified as the 



creation and annihilation operators, respectively, and satisfy the following commutation 
relations: 



\ak,a>i\ 



[4, af] = 0. K «!) = 5(3) (£ - 0- (3-10) 
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Using these relations, the equal-time two-point function for q is given as 
G(t,x;t,y) = (0\q{t,x)q{t,y)\0) 



71 



J ^^(^/^^fe^-^-w), (3.11) 



d /up 

AH inf J ~fay 

where physical momentum fc p hy S is related to the comoving one k as 

k phys = (a /a)k, (3.12) 

and Xphys = (a/ao)x. During the inflation, the present horizon scale is far outside of the 
horizon, and hence we are interested in the behavior of G(t,x;t,y) with |x p h ys — 2/ph ys | 3> 
Hfo f . In other words, we only need an information of the integrand with /c p h ys -C -f/inf- 
Denoting 

G(t,x;t,y) = f ^P^|g(A;)| 2 e^^-^, (3.13) 



q(k) for small k is given by 



m = 7 " . = (3.14) 

1 [—u + 1) smz/7r \2Hi n { J 2tt 

For the case where m 9 ^if inf , | (A: phys /2 J fZ' iri f ) - ^ -3 / 2 ^ | <C 1 since we are interested in modes 
with fcph y s *C H in f. (Notice that, in most of the inflation models, the COBE scale corre- 
sponds to ln(if in f//c phys ) ~ 50 — 60.) Thus, for the quintessence models with relatively large 
(effective) mass during the inflation, the primordial fluctuation in q is extremely suppressed 
relative to H in f. If the mass of the quintessence is light, however, v becomes close to 3/2 
and the suppression factor may become negligible. In this limit, q(k) is given by 



and hence q ~ H in f/27T for m q <C H m {. In order to avoid a significant suppression for 
the primordial perturbation, m q /H in{ ^ 0.1 — 0.2 is enough. In addition, notice that, for 
m q /H ini <C 1, the resultant q is expected to be almost scale independent; the only source 
of the scale dependence is a minor variation of the expansion rate H in f during the inflation. 
In the following discussion, we neglect the scale dependence of q. 

It is convenient to consider the ratio of the primordial value of q to that of the gauge- 
invariant variable ^ at the radiation-dominated universe 

r ^ih- (3 - 16) 

# 4 In our analysis, we only consider the case with scale-independent primordial fluctuation, and r q is 
treated as a scale-invariant quantity. When r q has a scale dependence via ^Sf and/or q, r q {k) for the present 
horizon scale becomes the most important since the quadrupole anisotropy is most strongly affected by 
the isocurvaturc mode as will be shown below. 
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where M* ~ 2.4 x 10 18 GeV is the reduced Planck scale, and \I> is defined as {^(x)^(y)) = 
Jrflnfc|^(A;)|V^-^. 

Solving the dynamics of the inflaton during the inflation, we obtain \l/ at the radiation- 
dominated era [|13|] 

where x is the inflaton field. Using the slow- roll condition for x, and also using q ~ H in f/27i, 
we obtain 

where V^ n f is the inflaton potential and V{ ni is its derivative with respect to the inflaton 
field. Q For example, for the chaotic inflation with Vi n f oc x p with p being an integer, r q is 
given by 

I 9PM * /o 1Q N 

4x(^cobe) 

where x(^cobe) represents the inflaton amplitude at the time when the COBE scale crosses 
the horizon. Numerically, we found r q ~ 0.3 — 0.6 for p = 2 — 10. Since the ratio r q 
generically depends on the model of the inflation, we treat free parameter in our 

analysis. 

Notice that, in our analysis, we treat \& and the fluctuation in the quintessence ampli- 
tude as random Gaussian variables. In particular, no correlation is assumed among these 
quantities. Thus, r q given in Eq. (|3.16|) should be understood as the ratio of the expecta- 
tion values of q and and we add adiabatic and isocurvature contributions in quadrature 
in calculating the CMB anisotropy. 

Now, we consider the evolution of the isocurvature fluctuation after the inflation. The 
evolution of the isocurvature mode is well described by the equation Eq. ( |3.2|) with ne- 
glecting the right-hand side (i.e., h = 0). This is because, when the energy density of the 
quintessence is negligible, the metric perturbation h is insensitive to the fluctuation in the 
quintessence and we can neglect the right-hand side of Eq. Q3.2p in studying the isocurva- 
ture perturbation. In the early universe, this is the case for the cosine-type quintessence. 
Let us first consider the case with large wavelength (i.e., /c p h ys "C H). In this case, the third 
term in Eq. ( |3.2| ) is irrelevant. By neglecting the right-hand side of Eq. ( |3.2|) , q behaves 
like a harmonic oscillator with effective mass-squared V"(Q) in the expanding background. 
Then, we can imagine two typical cases. If the effective mass is much smaller than the 

# 5 If Jq during the inflation, which is denoted as fg Di \ is smaller than the present value, r q is enhanced 

by the factor f Q /f^ ni) . This may happen if the quintessence is a pseudo Nambu-Goldstone boson; in that 
case /q is given by a vacuum expectation value of a (real) scalar field, which varies if the scalar potential 
is deformed during the inflation. 
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expansion rate of the universe, the slow-roll condition is satisfied for q. In this case, q 
keeps its initial value until the expansion rate of the universe becomes comparable to the 
effective mass of the quintessence field. This is the case for the cosine-type quintessence 
models in the early universe. 

As the universe expands after the inflation, the horizon expands more rapidly than the 
physical scale, and hence all the momentum scales relevant to our discussion eventually 
enter the horizon. Once the scale enters the horizon, the spatial derivative term cannot be 
neglected. The energy density of such modes behaves like that of relativistic matter, and 
hence it decreases as a -4 . Thus, once /c p h ys ~ H is realized, q(k) decreases as a -1 . 

So far, we have seen that the perturbation in the quintessence amplitude q can be gen- 
erated by several sources. If q is non-vanishing, the energy density of the quintessence field 
fluctuates. The energy density, pressure, and momentum perturbations of the quintessence 
are given by 

5 PQ = Qq + V'q, 5 PQ = Qq-V>q, (p Q + p Q ){v Q )i = - (^) Q^r, (3.20) 

respectively. These perturbations will affect the density perturbation of the universe and 
may change the behavior of the CMB anisotropy, as will be discussed in the following 
section. 

4 Numerical Results 
4.1 Outline of the Analysis 

To study the CMB anisotropy in models with quintessence, we solve the perturbed Einstein 
equation coupled to the equation of motion of the quintessence field (as well as evolution 
equations of other components like CDM, photon, baryon, and neutrinos). Then, we 
compare the theoretical prediction of the CMB anisotropy with observations. 

As a first step, the total density and pressure perturbations are derived by adding the 
quintessence contributions. Notice that these quantities affect the evolution of the metric 
perturbations. Then, we solve the evolution equations for the perturbations including the 
quintessence contributions. Evolution of Q(t) is simultaneously solved since it affects the 
expansion rate and the equation of state at each epoch. In particular, for the perturbed 
part of the equations, we modify the CMBFAST package ]14|] to include the quintessence 
contributions. Then, we calculate the CMB anisotropy for the Z-th multipole Ci, which is 
defined as 

(AT(x,7)AT(f,f)> 4H a + W-f). (4- 1 ) 

where AT(f, 7) is the temperature fluctuation of the CMB pointing to the direction 7, and 
the average is over the position x. 
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Since we have not specified any particular model of inflation which determines the 
normalization of the CMB anisotropy, over-all normalization of G\ is undetermined at this 
stage. We treat the normalization of the CMB anisotropy as a free parameter to make our 
study model- independent. Thus, in our analysis, the normalization of the CMB anisotropy 
is fixed so that the theoretical prediction has the best fit to the observational data. We 
denote the actual CMB anisotropy C\ as 



a = Nc h 



(4.2) 



where C\ is the normalization-free CMB anisotropy and N is the normalization factor which 
minimizes the \ 2 defined below. 

Since accurate measurements of the CMB anisotropy have been performed by COBE 
l| and MAXIMA [[|, we can constrain the quintessence models by 
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BOOMERanG 

comparing the theoretical prediction on C\ with observations. For this purpose, we calculate 
the goodness-of-fit parameter \ 2 = — 21nL, where L is the likelihood function. (We call 
this parameter as \ 2 i since it reduces to the usual x 2 -variable for a Gaussian.) In our 
statistical analysis, following Ref. [|1^], we use the offset lognormal approximation to derive 

X 2 as 



X 



BB' 



(4.3) 



where the summation is over the band powers obtained by COBE, BOOMERanG and 
MAXIMA. The quantity Z B hs contains informations from the observations as 



7obs 



hi(D B + x B ) 



(4.4) 



where D B is the observed band power in 5-th band and xb is the offset correction, and 
M BB , is given by 



M 



BB' 



M° B ,(D B + x B )(D B , + x B > 



where M BB , is the weight matrix for the band powers D B . (In Eq. 
over B nor B' is implied.) In addition, Z B h is written as 



(4.5) 

no summation 



7 til 



In 



Nj2 u cfBiQ + x B 



(4.6) 



where fsi is the filter function and u a is the calibration parameter.^] In our analysis, 
we neglect the calibration uncertainty in COBE and assume flat (top-hat) distributions of 
^BOOMERanG = 1 ± 0.2 and mmaxima = 1 ± 0.08 |TJ. For a given set of C h we vary N, 
■"BOOMERanG, and mmaxima to minimize the \ 2 parameter. In our numerical calculations, 
we use RADPACK package |18|] to calculate \ 2 i which is based on 24, 12, and 10 band powers 



from COBE, BOOMERanG, and MAXIMA, respectively. 

# 6 In our analysis, we assume that the calibration uncertainty u a has a flat distribution. Therefore, % 2 
given in Eq. (4.S) does not contain contributions from the calibration uncertainties. 
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Figure 4: The CMB angular power spectrum 1(1 + 1)Ci/2tt in models with cosine-type 
quintessence. The parameters we take here, /q = 1.8 x 10 18 GeV, A = 2.4 x 10~ 12 GeV 
(dashed line), A = 3.0 x 10~ 12 GeV (dotted line), A = 4.0 x 10~ 12 GeV (dash-dot line), 
A = 5.0 x 10~ 12 GeV (dash-dot-dot line). For comparison, we also show the cosmological 
constant case (solid line). Here, the cosmological parameters are taken to be h — 0.65, 
Q m = 0.3, VL^h 2 = 0.019, and the initial spectral index is n = 1. We also show the data 
points from COBE, BOOMERanG, and MAXIMA. (For COBE, we use the reduced data 



set given in Ref. [19]. 



4.2 Adiabatic Fluctuation 

Now, we are at the position to study the CMB anisotropy in models with quintessence. In 
this subsection, we consider the effects of the adiabatic perturbation so we take r q = 0. 

In Fig. [|, we first show the CMB angular power spectrum. Here, the cosmological 
parameters are taken to be h = 0.65, Q m = 0.3, Q^h 2 = 0.019, and the initial spectral 
index is n = 1. We assume that there are no tensor mode contributions. The normalization 
factor is chosen such that x 2 is minimized. 

Some of the interesting features of the CMB angular power spectrum are discussed in 
order. First, let us consider the locations of the acoustic peaks. The locations of the peaks 
depend on two quantities, the sound horizon at last scattering and the angular diameter 
distance to the last scattering surface. Approximately, the Z-th multipole picks up scales 
around / ~ kreij^) where r^r*) is the angular diameter distance to the last scattering 



surface |[20|| . The n-th peak in the temperature power spectrum is located at the scale k n 
which satisfies k n r s (T*) = nir, where r s (rA is the sound horizon at last scattering. So the 
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location of n-th peak in the I space is estimated as 



L ^ — 7— v-tot. (4.7) 

Since the behavior of the cosine-type quintessence is almost the same as that of the cos- 
mological constant until very recently, the sound horizon at last scattering is the same in 
both cases. The angular diameter distance in the quintessence model is, however, different 
from that in the ACDM models. Since the quintessence models provides larger total energy 
density of the universe than the ACDM models in the earlier epoch, the angular diame- 
ter distance in the quintessence model becomes smaller than that in the ACDM models. 
As we can see from Fig. |], the location of the peaks is shifted to lower multipole I for 
the quintessence models. If we take a parameter in the oscillatory region, this feature 
becomes more prominent. As A becomes larger, more energy density exists in the early 
universe since some fraction of the energy density of the quintessence damps away during 
the oscillation. 

Next let us consider the height of the acoustic peaks. Since the energy density of 
the quintessence becomes dominant when z ~ 0(1), the late time integrated Sachs- Wolfe 
(ISW) effect enhances low multipoles. Such an enhancement may be more effective in the 
quintessence models than in the ACDM models since, in the quintessence case, the "dark 
energy" (i.e., the quintessence) may dominate the universe earlier than in the ACDM case. 
As a result, with the COBE normalization, the height of the first peak becomes lower. On 
the contrary, since the quintessence becomes the dominant component of the universe only 
at later epoch, pattern of the acoustic oscillation before the recombination does not change 
compared to ACDM models. Therefore, ratios of the height of the first peak to those of 
higher peaks are the same as ACDM models. 

There is another point which should be addressed. If we take a parameter in the deep 
oscillatory region, we can see a bump in the low multipole region. At low multipoles, the 
late time ISW effect is important for the temperature fluctuations in ACDM models or in 
the quintessence models. The ISW effect is originated from the decay of the gravitational 
potential at the time of the quintessence (or the cosmological constant) domination. Before 
the quintessence field starts to oscillate, the equation of state of the quintessence field 
is almost —1, and hence this epoch is like the cosmological-constant-dominated epoch. 
After the quintessence starts to oscillate, however, ujq approaches to and hence the 
quintessence field behaves like a matter component. Then, the gravitational potential \l/ 
takes a constant value again, and the ISW effect becomes ineffective. Consequently, the 
bump in low multipoles shows up. Notice that the location of this bump corresponds to 
the transition scale from the slow-roll epoch to the oscillatory epoch. 

Now we discuss constraints on the cosine-type quintessence models from the obser- 
vations of COBE, BOOMERanG and MAXIMA. Following the prescription given in the 
previous subsection, we calculate the goodness-of-fit parameter x 2 as a function of Jq and 
A. Based on this x 2 variable, constraints on the parameter fq and A are shown in Fig. 
H Here, we show the region consistent with x 2 < 60.5 and x 2 — 68.7, which correspond 
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Q Q <0.7 



0.5 1 1.5 2 2.5 3 
f Q ( 10 18 GeV) 

Figure 5: Constraints on the parameters A and /q. The lightly shaded regions are for 
X 2 < 68.7, and the darkly shaded region is for \ 2 < 60.5. The cosmological parameters are 
taken to be h = 0.65, Q m = 0.3, VL^h 2 = 0.019, and the initial spectral index is n — 1. The 
Points A — E will be used as representative points in the later discussion. 

to 95 % and 99 % C.L. allowed region for the x 2 -statistics with 44 degrees of freedom, 
respectively. 

When we take a parameter in the oscillatory regions, the quintessence field becomes 
dominant component of the universe at earlier epoch. Namely the late time ISW effect 
becomes large and the angular diameter distance to the last scattering surface becomes 
smaller as mentioned before. Therefore, if A is significantly large, the late time ISW effect 
enhances angular power spectrum at low multipoles. As a result, the heights of the acoustic 
peaks are suppressed relative to C\ with small I. 

4.3 Isocurvature Mode 

Now, let us consider the effect of the isocurvature fluctuation which may be generated in 
the very early universe (like during the inflation). 

Before discussing the CMB anisotropy generated by the isocurvature mode, let us first 
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k(h Mpc" 1 ) 

Figure 6: q as a function of k. q is normalized to unity at large scale (small k). Here, the 
model parameters are taken to be J'q = 5.0 x 10 17 GeV and A = 2.3 x 10~ 12 GeV (solid 
line), f Q = 5.0 x 10 17 GeV and A = 3.5 x lO" 12 GeV (dashed line) . 



study the behavior of the fluctuation in the quintessence amplitude q. Since the isocur- 
vature contribution to C\ does not interfere with the adiabatic one, studying the purely 
isocurvature case is enough to understand the effect of the isocurvature perturbation. 

In Fig. H, we show the scale dependence of q. As one can see, q{k) changes its behavior 
at the scale which is slightly smaller than the present horizon scale: q(k) is almost constant 
for small k while q(k) is suppressed for large k. This behavior is understood as follows. 
When the slow-roll condition is satisfied for the quintessence, fluctuation for the scale larger 
than the horizon scale takes its initial value. On the contrary, for A; p h y s ~ H, q behaves 
like an oscillator with frequency ~ fcphys, and the fluctuation red-shifts as the universe 
expands. Due to the red shift, the magnitude of q is approximately proportional to k~ 2 
and for modes entering the horizon in the matter-dominated and radiation-dominated 
universe, respectively. If the slow-roll condition is satisfied until the present universe, q{k) 
for fcphys ~ Hq is suppressed. 

If m e ff ~ A 2 / /q ^ H , however, the quintessence field may start to oscillate at some 
stage of the universe. (We call the horizon scale at this epoch as 1/k^l.) Once the 
quintessence starts to oscillate with frequency m e g, then all the modes universally damps 
as far as A; p h y s ~ rn c g- Thus, for fc p h ys ~ k^l , q is independent of k. On the contrary, 
modes with /c p h ys ~ kphyl acquire extra suppression as we discussed, and q for such a small 
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Figure 7: 5q as a function of k. 5q is normalized to unity at large scale (small k). Cosmo- 
logical and model parameters are the same as Fig. ^|. 

scale becomes smaller. In our analysis we only consider cases where 1 / k^l is close to the 
present horizon scale. In those cases, the isocurvature fluctuation in the quintessence field 
dominantly affects C/ with small I, as we will see below. 

In Fig. 0, we also plot the scale dependence of the energy density perturbation in the 
quintessence field normalized by pg 



5_pg = Q'q + V'q 
PQ PQ 



For the isocurvature contribution, the energy density perturbation is dominated by that of 
the quintessence and hence it induces the metric perturbation ^. This becomes a source 
of the CMB anisotropy through the Sachs- Wolfe effect. When the slow- roll condition is 
satisfied for the quintessence, 5q(/c) is approximated as 



5, 



Q 



YL 

PQ 



2q 



9H 



+ q 



^phys 



(4.9) 

A 2 //q, as a consequence 



and it changes its behavior at the scale /c p h ys ~ Ho or at k v 
of the scale dependence of q. 

Now, we consider the CMB anisotropy in the case with the isocurvature mode. We 
calculate the CMB anisotropy for various cases, and in Table [1], we show the quadrupole 
C 2 normalized by C 10 . 
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T'q = r q = 0.5 r q = 1 r g = 1.5 r q = 2 



A 


1.31 


1.31 


1.32 


1.33 


1.34 


B 


1.45 


1.46 


1.48 


1.51 


1.56 


C 


1.15 


1.26 


1.61 


2.18 


2.98 


D 


0.84 


1.33 


2.80 


5.25 


8.69 


E 


0.62 


1.93 


5.86 


12.41 


21.58 



Table 1: C 2 /C w for several values of r q . We take f Q = 5x 10 17 GeV, and (A) A = 2.3 x 10~ 12 
GeV, (B) A = 2.5 x 10~ 12 GeV, (C) A = 3.0 x 10" 12 GeV, (D) A = 3.5 x 10~ 12 GeV, and 
(E) A = 4.0 x 10~ 12 GeV. Notice that the Points A - E are indicated in Fig. g. 



If we limit ourselves to the parameter region which is consistent with the COBE, 
BOOMERanG, and MAXIMA observations with simple scale-invariant primordial fluc- 
tuation, effect of the isocurvature mode is quite small as far as r q ^z 1. This is because, 
in such cases, there is a severe upper bound on A to suppress the late time ISW effect 
which enhances Ci with small I. As a result, the quintessence field cannot dominate the 
universe when z ^> 1. Then, the isocurvature fluctuation in the quintessence density also 
becomes a minor effect until very recently. As one can see in Table [j], for the best-fit 
value of A (i.e., for the Point A given in Fig. |5|), the enhancement of C 2 is about 2 % 
even for tq = 2. If we consider larger value of A, effect on C 2 is more enhanced. For 
(/ Q ,A) = (5 x 10 17 GeV, 3.0 x 10~ 12 GeV) (i.e., for the Point C given in Fig. |, which is 
allowed at 99 % C.L.), we calculate the CMB angular power spectrum and the result is 
given in Fig. |8|. In this case, C 2 can be enhanced by the factor 2.6 if tq = 2. 

One should note that angular power spectrum of the CMB anisotropy strongly depends 
on the primordial spectrum of the fluctuations. Thus, the constraint on the J'q vs. A 
plane is sensitive to the scale-dependence of the primordial adiabatic fluctuation which is 
determined by the model of the inflation. Therefore, if we adopt a possibility of a non- 
trivial scale dependence of the primordial fluctuation, the constraint on the fq vs. A plane 
given in the previous section may be relaxed or modified. If this is the case, larger value of 
A may be allowed and the energy density of the quintessence field may become significant 
at earlier stage of the universe. Thus, we also consider such cases. In Fig. || we plot Ci 
normalized by C 10 for (/ Q ,A) = (5 x 10 17 GeV, 3.5 x 10~ 12 GeV). In addition, in Table 
|l], the ratio C 2 /Ciq is shown for the Points D and E. As one can see, effect on C 2 is much 
more significant than the previous case. 

At low multipoles the uncertainty of the CMB data is dominated by "cosmic variance." 
For C 2 which is mostly enhanced by the isocurvature fluctuations, the cosmic variance 



gives y2/5 ~ 60 % error to observation.^ 7 ! Thus, it is difficult to see the effect of the 

V 



isocurvature mode when r q ^ 1 if the quintessence field does not oscillate so much until the 



# 7 The measurement of the CMB polarization toward many clusters would allow some further reduction 



of the cosmic variance |21 
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Figure 8: Q/C 10 for r q = 0, 0.5, 1, 1.5, 2, 3, and 4 from below with f Q = 5x 10 17 GeV 
and A = 3.0 x 10~ 12 GeV. 
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Figure 9: Ci/C w for r g 
A = 3.5 x 10- 12 GeV. 



0, 0.5, 1, 1.5, and 2 from below with f Q = 5 x 10 17 GeV and 
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present epoch. However, if a large value of A is possible, the isocurvature effects may be 
detectable even with r q ^ 0.5, which is realized, for example, in the chaotic inflation model 
with Vi n f oc x p with 8 (see, for example, Point E in Table []]). 



5 Conclusions and Discussion 

We have studied the CMB anisotropies produced by cosine-type quintessence models. In 
particular, effects of the adiabatic and isocurvature fluctuations have been both discussed. 

For purely adiabatic fluctuations with scale invariant spectrum, the existence of the 
quintessence suppresses the relative height of the first acoustic peak of the angular power 
spectrum compared with the ACDM case. This is because, in the quintessence models, 
the "dark energy" due to the quintessence may dominate the universe earlier than the 
cosmological constant case, and hence the late time ISW effect becomes more effective. As 
a result, the CMB anisotropy for large angular scale is more enhanced, which relatively 
suppresses the height of the acoustic peaks. Because of this effect, we have seen that 
the CMB data from COBE, BOOMERanG and MAXIMA have imposed the stringent 
constraint on the model parameters of the quintessence models. We have also seen that 
the location of the the first acoustic peak shifts to lower multipole / compared with ACDM 
models. 

In the case of the cosine-type quintessence models, the quintessence field has a negligible 
effective mass during inflation and hence its amplitude may acquire sizable fluctuation 
due to the quantum fluctuation in the de Sitter background. Such a fluctuation becomes 
isocurvature fluctuation. We have shown that the isocurvature fluctuations have significant 
effects on the CMB angular power spectrum at low multipoles in some parameter space, 
which may be detectable in future satellite experiments. This signal may be used to test 
the cosine-type tracker model, combining with the global shape of the CMB angular power 
spectrum. In particular, enhancements at low multipoles do not exist in the minimal 
ACDM models. 

In the tracker-type models, the CMB angular power spectrum may be also affected by 
the isocurvature mode f22fl . Contrary to the cosine-type case, however, tracker field starts 



to evolve in the early universe and hence its effective mass is as large as the expansion rate 
of the universe. As a result, fluctuation in the tracker field damps while the tracker field 
follows the attractor solution. In addition, in the tracker case, the flatness of the tracker 
potential is not guaranteed by any symmetry during the inflation. (Notice that for the 
cosine- type case some symmetry may keep the flatness of the potential.) As a result, in the 
de Sitter background, the tracker field may acquire effective mass as large as H ini during 
the inflation. This is the case for, for example, quintessence models based on supergravity 
with minimal Kahler potential. If such a large mass exists, fluctuation generated during 
the inflation is damped, as can be seen in Eq. (|3.15|) . These effects may drastically suppress 



the signal of the isocurvature mode. Detailed analysis of the tracker models will be given 



in elsewhere 23 . 
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Since very accurate measurements of the CMB anisotropy are expected in the near 



future, in particular by the MAP |24| and PLANCK |2fJ experiments, we will have more 
stringent constraints on the quintessence models as well as on the ACDM models. Impor- 
tantly, the quintessence and ACDM models may have different predictions on the shape of 
the CMB angular power spectrum, and some of the models may be confirmed or excluded 
once better observations of the CMB anisotropy become available. In particular, in the 
cosine-type quintessence models, we have seen that the isocurvature perturbation may en- 
hance the CMB angular power spectrum at low multipoles (in particular, C2), which may 
be an interesting signal of the cosine-type quintessence models. 
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